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Abstract: Numerical investigations are presented for a multilayer frequency selective surface with Koch 
fractal (levels I and 2) conducting patch elements. The structure investigated is obtained using two FSS 
screens separated by an air gap layer. For the proposed investigation were used three different values an 
air gap height. The results obtained using the numerical method were compared with other technique and 
using the commercial software Ansoft Designer™. A good agreement was observed in terms of the 
bandwidth. 
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I. Introduction 

A periodic surface can be defined basically as a set of identical elements arranged in two or three 
dimensions forming an infinite array. A periodic array formed by conductive patch elements or opening 
elements is known as Frequency Selective Surface (FSS). FSS are periodic structures that can provide 
frequency filtering to incoming electromagnetic waves, and their frequency response is entirely determined by 
the geometry of the structure in one period called a unit cell [1,2]. 

Traditional FSS has a long history of development and have been investigated over the years for a 
variety applications, e.g., frequency filters or diplexers in high performance reflector antenna systems, 
advanced radome designs, and smart surfaces for stealth applications [3]. 

The question of operating bandwidth is one of important problems in FSS theory, in some applications 
a multiband frequency response is desired; however in other applications a wideband frequency response is 
more preferred. FSS structures have been successfully proven as a mean to increase the communication 
capabilities of satellite applications. In space mission (Voyager, Galileo and Cassini), for example, the use of 
dual antenna system with FSS reflectors has made it possible to share the main reflector among different 
frequency bands, therefore FSS with multiband responses have been studied by several researchers. But, in 
some applications, such as separating successive frequency bands, FSS are required to have a wide 
transmission or reflection bands and rapid rolloffs [4-10]. 

Experimental and numerical investigations on cascading or multilayer FSS have been done, using full 
and simple (approximate methods) [11-14]. This paper presents a numerical investigation of multilayer or 
cascaded FSS using a simple but efficient method used to predict the frequency response (transmission 
coefficient) for these structures. In this work was analyzed a multilayer FSS structure, using conducting Koch 
fractal elements. This structure was first proposed in [15] and consists of two FSS, called structure 1 and 
structure 2 mounted on a dielectric isotropic layer, separated by an air gap layer. The dielectric substrate used 
was the RT-Duroid 3010, with 1.27 mm of height and relative permittivity equal to 10.2 and unit cell 
periodicity (7^ = T y ) equal to 10 mm. We used the commercial software Ansoft Designer™ to obtain the 
individual scattering parameters of each one FSS structure. To validate the results a comparison with results 
from other numerical method is performed. 

II. Fss Theory 

Through the use of the spectral domain approach to analyzing the response of an FSS, which is 
assumed to have infinitesimal thickness, one can determine the electromagnetic fields on a plane z = z m see 
Fig. 1 and Fig. 2, These electric fields are expressed in terms of a discrete spectrum of plane waves known as 
Floquet harmonics. 
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Fig. 1: Front (superior) view of a FSS. 
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Fig. 2: FSS Cascading, side view. 

Cwik in [16] has defined a set of scattering parameters that incorporates the vector nature of electromagnetic 
fields as follows [1]: 



S n (m,nJJ) 
S 12 (m,nJJ)-- 
S 21 (m 9 n,iJ) = 
S 22 (m,nJJ)-- 



<Jp (ij) 



slP inc (i,j) 

where zr is the above interface and zl is the below interface in Fig. 1, ie, zr = 0 and zl = z p+q , are Floquet 
voltage waves given by [1]: 

V (±,±) (Km > k yn ' Ki > k yj >Z R )=f(m, It, 1, j, Z) (±±) ^P(m,n), 
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With, 

P(m,n) = (kl+k 2 yn )Y mn , 

Where: 

y Y mn 

mn . ' 

JCQ/U 

The terms K xm and K yn (wave numbers of Floquet special harmonics) are defined as: 

K m = Y m+k '^ 

X 

K n =Y n + k 7> 

y 

The incident harmonics are functions of elevation angle and the azimuthal angle, which are given as [1]: 

K c =Ksen6 inc cos4 c , 

k™ =k o sen0 inc sen0 inc , 

Referring to the coordinate system of Fig. 1, can be seen that the notation (+, +) associated with the 
definition for the voltage waves indicates the direction, both +z or -z, of the incident and scattered energy in a 
FSS system. Regarding the scattering parameters, it is more convenient to express the potential vector, used in 
the definition of the Floquet voltage waves in terms of the electric fields transform. 
The potential vector, expressed in terms of total electric fields, can be written as [1]: 



i(k E ± -k E ± ) 

f(m,n,i,j,z) (±±) = - yn :r 7 y - } , 



k z + k z 

K xm ^ K yn 

The total electric field can be found as follows [1]: 

E; , (k ,k ) = E: , +E ( , 8(k -k )S(k -k ), 

E7 , (k ,k ) = E7 , +E ( , S(k -k )S(k -k ), 

( x >y) t otai y xm ' y n ' ( x >y)s ( x >y)tmn S v Xl xm ' v yj yn'' 

Where E (re f t trans) are the reflected and transmitted electric fields, calculated in z = 0 and z = z p+q . We can write 
the scattering parameters in terms of the electric field evaluated at the reference planes above (zr) and below 
(zl) of a FSS: 



S n (m,n,i,j) = 



' Y u pl+k'M+k 2 ) 



S l2 (m,n,i,j)=l-^ 



Y j(k E~ -k E~ ) 

1 ntn Jy yn x, oM xm y total ' 



S 2l (m,n,i,j) = 



WZ ^KKt,-K m KJ 



IY j(k E~ -k E + ) 



Y ij pL+k 2 yn )(k 2 xi +k 2 yj ) 
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III. Numerical Method 

Since scattering parameters (scattering matrix) of finite dimension for a FSS are computed, many 
analytical procedures are available for obtaining multilayer composite representation. It is possible to make 
direct use of the scattering matrix along with the following equation to obtain the representation of a system 
with two FSS (multilayer or cascaded FSS), each FSS (individual structure) is viewed as a subsystem [1]. The 
scattering matrix for the multilayer FSS, using two FSS, is given by [17]: 



s c = 



(1) 



+ S$PTS^PS 21 

A2) D/ C(l) , C(l) DTC( 2 ) DC(1) > 



c(l) 

°11 



(2) 



S£PTS U 



(2) 



22 12 



+ S. 



(2) 
22 



Where , , and are the scattering parameters which represent the first subsystem, namely, 



the first FSS, , S$ , and are the scattering parameters associated with the second subsystem, 
and 5 C is the scattering matrix for the system composed by two FSS (multilayer FSS). 
And Tis given by [17]: 

T = (i-s£ ) ps%py\ 



*(2) 



: (2) 



P is a diagonal matrix that has e ( ~ j2kd) as its elements, k is the wave number [45], in this work k = k 0 (because 
the multilayer FSS is formed by two single FSS separated by an air gap layer), and / is an identity matrix. 
To validate the results, that will be presented, a comparison with results from other numerical method called 
"One Mode Interaction Technique" that can be seen in [17] and obtained by Method of Moments (MoM) 
using the commercial software Ansoft Designer™ is performed. The numerical results were computed using 
the commercial software MATLAB™ (MATrix LABoratory). 

Using the One Mode Interaction (OMI) technique, the transmission (C T ) and reflection (C R ) coefficients for the 
multilayer FSS, in this case formed by two single FSS, are [18]: 



TT 

1 l 1 2 



(-j2kd 1 ) ' 



C R =R l + 



S-j2kd x ) 



\-R x R 2 e 



(-j2kd l ) 



Where dj is the spacing between the two structures. T h Rj are the transmission and reflection coefficients for 
the first FSS, and T 2 , R2 are the transmission and reflection coefficients for the second FSS. 

IV. MULTILAYER FSS STRUCTURE 

The multilayer FSS used in this work consists of two FSS screens called structure 1 and structure 2, 
respectively, each one using Koch fractal patch elements printed on a dielectric substrate separated by an air 
gap layer. The first structure with Koch fractal level 1 is defined as structure 1 and the second FSS screen with 
Koch fractal level 2 is defined as structure 2, as can be seen in Table 1 with their respective resonant 
frequencies if) and bandwidths (BW). The process to generate the Koch fractal elements used in this work can 
be seen in [14]. The element shapes and the multilayer structure considered in our investigation are shown in 
Fig. 3 and Fig. 4, respectively. 



Table 1 - Identification of the elements 



Structure 


Element type 


/(GHz) 


BW(GHz) 


1 


Koch fractal level 
1 


9.115 


1.951 


2 


Koch fractal level 
2 


8.333 


1.618 
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Fig. 3: Koch fractal shape elements. 
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Fig. 4: Multilayer Structure. 



V. Results 

In order to validate the results obtained with the numerical method presented before, they are 
compared with results obtained with the OMI technique and the results obtained using the commercial 
software Ansoft Designer™, the latter uses the full wave method, MoM, to compute its results. 

Fig. 5 illustrates the results obtained with the numerical method for the case with an air gap height 
equal to 2.0 mm. For the multilayer FSS with an air gap height equal to 2.0 mm, the results obtained using the 
numerical method presents a bandwidth, for 20 dB insertion loss reference, equals to 3.420 GHz, a good 
agreement is observed between the results. 

For the multilayer FSS with an air gap height equal to 4.0 mm, the results obtained using the 
numerical method presents a bandwidth, for 20 dB insertion loss reference, equal to 3.360 GHz, 
approximately. Fig 6. illustrates a comparison between the results. 

Using an air gap height equal to 6.0mm, the numerical method presents a bandwidth, for the 20 dB insertion 
loss reference, equal to 3.290 GHz, approximately. 
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Fig. 5: Transmission coefficient for the multilayer FSS structure with an air gap height equal to 2.0 mm. 
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Fig. 6: Transmission coefficient for the multilayer FSS structure with an air gap height equal to 4.0 mm. 
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Fig. 7: Transmission coefficient for the multilayer FSS structure with an air gap height equal to 6.0 mm. 

VI. Conclusion 

In this paper, a multilayer FSS was investigated using a simple numerical method. This multilayer 
structure was formed by two FSS screens separated by an air gap layer, each one FSS screen using a 
conducting patch element with fractal geometry in the unit cell. The results were obtained for the numerical 
method and compared with other numerical technique called "One Mode Interaction". Moreover, both were 
compared with the results obtained using the commercial software Ansoft Designer™ for different values of 
spacing between the FSS screens. The numerical method is efficient and can be used in conjunction with other 
methods, like full wave methods such like: Finite Element Method, Wave Concept Iterative Procedure, and 
Equivalent Circuit Model, among others. For all the cases considered, the frequency responses have 
characteristics of a high-pass filter with a very large rejection band. A good agreement between the results was 
observed in terms of bandwidth. 
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